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We present a detailed microscopic theory of the conserved spin current which is introduced by us 
[Phys. Rev. Lett. 96, 196602 (2006)] and satisfies the spin continuity equation even for spin-orbit 
coupled systems. The spin transport coefficients o~i v as a response to the electric field are shown to 
consist of two parts, i.e., the conventional part affl and the spin torque dipole correction o s J v . As one 
key result, an Onsager relation between erf v and other kinds of transport coefficients are shown. The 
expression for in terms of single-particle Bloch states are derived, by use of which we study the 
conserved spin Hall conductivity in the two dimensional hole gas modeled by a combined Luttinger 
and SIA Rashba spin-orbit coupling. It is shown that the two components in spin Hall conductivity 
usually have the opposite contributions. While in the absence of Rashba spin splitting, the spin 
Hall transport is dominated by the conventional contribution, the presence of Rashba spin splitting 
stirs up a large enhancement of the spin torque dipole correction, leading to an overall sign change 
for the total spin Hall conductivity. Furthermore, an approximate two-band calculation and the 
subsequent comparison with the exact four-band results are given, which reveals that the coupling 
between the heavy hole and light hole bands should be taken into account for strong Rashba spin 
splitting. 

PACS numbers: 73.20.-r, 71.15.-m 



I. INTRODUCTION 



Spintronics which combines the basic quantum me- 
chanics of coherent spin dynamics and technologi- 
cal applications in information processing and storage 
devicesi' 2 -' 3 -, has grown up to become a very active and 
promising field in condensed matter. One central issue 
in spintronics is on how to generate and manipulate spin 
current as well as to exploit its various effects in a vari- 
ety kinds of systems, ranging from ferromagnetic metals 
to semiconductor paramagnets. In the ideal situation 
where spin is a good quantum number, spin current is 
simply defined as the difference between the currents of 
electron carried by the two spin states. This concept of 
the spin current has served well in early studies of spin- 
dependent transport effects in metals. The ubiquitous 
presence of spin-orbit coupling inevitably makes the spin 
non-conserved, but this inconvenience is usually put off 
by focusing one's attention within the so-called spin re- 
laxation time. In recent years, it has been found that 
the extrinsic or intrinsic spin-orbit coupling can provide 
a route to generate transverse spin current in ferromag- 
netic metals^ or semiconductor paramagnets^^ by the 
driving of an electric field. The fundamental question of 
how to define the spin current properly in the general 
situation then needs to be answered. In most of previous 
studies of bulk spin transport, it has been conventional 
to define the spin current simply as a combined ther- 
modynamic and quantum-mechanical average over the 
symmetric product of spin and velocity operators. Un- 



fortunately, no viable measurement is known to be pos- 
sible for this spin current. The recent spin-accumulation 
experiments ^ 10 ' 11 do not directly determine it, and there 
is no deterministic relation between this spin current and 
the boundary spin accumulation. 

In fact, the conventional definition of spin current suf- 
fers critical flaws that prevent it from being relevant to 
spin transport and accumulation. First, this spin cur- 
rent is not conserved, rendering it useless in describing 
a true "current" ; Second, this spin current does not nec- 
essarily vanish in insulator s 12 ' 13 , and in thermodynamic 
equilibrium^, so it is disqualified as a true transport cur- 
rent corresponding to spin accumulation; Finally, there 
does not exist a mechanical or thermodynamic force in 
conjugation with this current, so it cannot be fitted into 
the standard near-equilibrium transport theory. The last 
issue in particular makes the direct measurement of the 
conventional spin current difficult if not impossible. For 
instance, because the conventional spin transport coeffi- 
cients cannot be associated with other transport coeffi- 
cients via Onsager relation^, they cannot be measured 
by linking to other transport phenomena. 

These issues were addressed in our last brief report^, 
where we have established a proper definition of spin cur- 
rent free from all the above difficulties. The new spin cur- 
rent is given by the time derivative of the spin displace- 
ment (product of spin and position observables), which 
differs from the conventional definition by a torque dipole 
term. The torque dipole term is first found in a semiclas- 
sical theory^ 7 -, whose impact on spin transport has been 
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further analyzed to assess the importance of the inverse 
spin Hall effect^, i.e., the charge Hall effect driven by 
a spin force. The new spin transport coefficients <r* 
are shown to consist of two parts, i.e., the conventional 
part a a ® v and the spin torque dipole correction cr s J v . As 
one key result, the Onsager relation between <j*„ v and 
other kinds of force-driven transport coefficients has been 
shown. Note that the other alternative definitions of spin 
current have also been proposed recentl y 18 i 19 ' 20 i 21 ' 22 . 



In this paper, a detailed quantum- mechanical linear re- 
sponse theory of this conserved spin current is addressed, 
which will be shown itself a necessary supplement as 
well as an enlightening illustration for our previous brief 
report—. In particular, a general Kubo formula for the 
spin transport coefficients a s ^ v in terms of single-particle 
Bloch states is given in this paper. Then we use our 
formula to study the conserved spin Hall conductivity 
in the two dimensional hole gas (2DHG) modeled by a 
combined bulk Luttinger and space-inversion asymmet- 
ric (SIA) Rashba spin-orbit coupling. It is shown that 
the two components in spin Hall conductivity usually 
have the opposite contributions. While in the absence 
of Rashba spin splitting, the spin Hall transport is dom- 
inated by the conventional contribution, the presence of 
Rashba spin splitting stirs up a large enhancement of the 
spin torque dipole correction, leading to an overall sign 
change for the total spin Hall conductivity. Furthermore, 
an approximate two-band calculation and the subsequent 
comparison with the exact four-band results are given, 
which reveals that the coupling between the heavy hole 
and light hole bands should be taken into account for 
strong Rashba spin splitting. 



Our paper is organized as follows. In Sec. II we give 
a brief review of the conserved spin current defined in 
Ref The intuitive arguments begin with a search for 
spin density continuity equation by taking into account 
the spin torque dipole correction. In Sec. Ill we clar- 
ify the necessity to include the spin torque dipole cor- 
rection into the total spin current by initiating a linear 
response analysis. That is, only when this spin torque 
dipole correction is included, can the proper Onsager re- 
lations between the spin transport coefficients and the 
other transport coefficients be established, which turn 
out to be essential to endow the spin current a driving 
force. In Sec. IV we present a Kubo formula for the 
conserved spin transport coefficients in terms of single- 
particle eigenstates. In Sec. V we give a detailed discus- 
sion of the spin Hall effect in the two dimensional hole 
gas by use of our general formulae. The analytic and 
numerical calculation are given both in weak and strong 
SIA spin-orbit coupling regimes and show very different 
features concerning the interplay between the two com- 
ponents in the total spin Hall conductivity. Finally, in 
Sec. VI we present our conclusions. 



II. SPIN CONTINUITY EQUATION AND 
INTRODUTION OF CONSERVED SPIN 
CURRENT 

The conventional definition of the spin current is the 
expectation value of the product of spin and velocity 
operators. In the case of spin-polarized magnetic sys- 
tems, this spin current is simply reduced to the difference 
between the spin-up and spin-down electrical currents. 
In the case of spin-orbit stronly coupled semiconductor 
systems, where time-reversal symmetry ensures no bulk 
spin polarization, unfortunately, no viable measurement 
is known to be related with the conventional spin current. 
This can be straightforwardly shown by writing down the 
continuity equation for spin density, 



dSz 
dt 



+ V • J s = T 2 



(1) 



The spin density for single-particle (spinor) state ip(r) 
is defined by S z (r) = ip'(r)s z ip(r), where s z is the spin 
operator for a particular component (z here, to be spe- 
cific). The spin current density is given by the conven- 
tional definition J s (r) = Rei/^ (r)^{v, s z }-0(r), where v 
is the velocity operator, and {,} denotes the anticom- 
mutator. The right hand side of the Eq. (TT]) is the spin 
torque density defined by T z (r) = Keip^ (r)r^(r), where 
f = ds z /dt = (l/ih)[s z , H], and H is the Hamiltonian 
of the system. The presence of the torque density T z re- 
flects the fact that spin is not conserved microscopically 
in systems with spin-orbit coupling. As a consequence, 
the only knowledge about the experimentally measured 
variation of spin density in space-time is not sufficient 
to determine the conventional spin current J s , and vice 
versa, due to the unique presence of the spin torque den- 
sity T z in spin-orbit coupled systems. 

One promising choice^ to remedy this oblique rela- 
tionship between spin density and spin current is to for- 
maly move the torque density term to the left hand side 
of Eq. (1) and absorb it in the divergence term. The 
physical reason is that, due to symmetry reasons, it of- 
ten happens that the average torque vanishes for the bulk 
of the system, i.e., (1/V) / dVT z (r) = 0. Thus one can 
write the torque density as a divergence of a torque dipole 
density, 

T 2 (r) = -V-P T (r). (2) 
Moving it to the left hand side of Eq. {T]), one has 
dS 



dt 



+ V-(J s +P r )=0, 



(3) 



which is in the form of the standard sourceless continuity 
equation. This shows that the spin is conserved on av- 
erage in such systems, and the corresponding transport 
spin current is: 



\T s — J s P*r ■ 



(4) 
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Note that there is still an arbitrariness in defining the 
effective spin current because Eq. Q does not uniquely 
determine the torque dipole density P r from the cor- 
responding torque density T z . However, this ambiguity 
can be eliminated by imposing the physical constraint 
that the torque dipole density is a material property that 
should vanish outside the sample. This implies in partic- 
ular that JdVP T = - JdVrV ■ P T = J dVrT z (r). It 
then follows that, upon bulk average, the effective spin 
current density can be written as J s = Kct/j* (r)J s tp(r), 
where 



Js = 



d(fs 2 
dt 



(5) 



is the effective spin current operator. Compared to the 
conventional spin current operator, it has an extra term 
r(ds z /d£), which accounts the contribution from the spin 
torque. 

The conservation of the new spin current allows one to 
consider spin transport in the bulk without the need of 
laboring explicitly a spin torque (dipole density) which 
may be generated by the electric field. Thus the spin 
transport in spin-orbit coupled systems can be treated in 
a unified way, whatever the spin-orbit coupling strength 
is weak or strong. For example, it has been customary to 
link spin density and spin current through the following 
phenomenological spin continuity equation, 

dS? S? ,„ N 

1^+V-.7 S = -, 6 

at t s 

where r s is the spin relaxation time, and the spin current 
has the form J s = oE — D S WS Z . This makes sense only 
if our new spin current is used in the calculation of spin 
conductivity a, otherwise an extra term of field-generated 
spin torque must be added. 

Equation ([6]) now can serve as the basis to determine 
the spin accumulation at a sample boundary. Consider a 
system having a smooth boundary produced by a slowly 
varying confining potential. We assume that the length 
scale of variation is much larger than the mean free path, 
so that the above continuity equation may be applied 
locally. By integrating from the interior to the outside 
of the sample boundary, we obtain a spin accumulation 
per area with S z = Jg t s . Thus one can see that 
for the generic class of smooth boundaries, there is a 
unique relationship between spin accumulation and the 
conserved spin current J s , instead of the conventional 
spin current J s . The other kinds of boundary conditions 
have also been discussed in previous wor k 16 i 23 ' 24 to clarify 
the relationship beween spin accumulation and bulk spin 
current. 



III. SPIN HALL CONDUCTIVITY AND ITS 
ONSAGER RELATION WITH INVERSE SPIN 
HALL EFFECT 

Defined as a time derivative of the spin displacement 
operatore fs z , the new spin current has a natural conju- 



gate force, i.e., the spin force, F s . To show this, one can 
consider the system exposed to an inhomogeneous Zee- 
man magnetic field along z-axis. The resulting perturba- 
tion can be modeled as V(r) = j*/iflB(r)s 2 with Bohr 
magneton /ib and effective magnetic factor g*. Suppose 
that the inhomogeneous Zeeman field is smoothly vary- 
ing in space around zero. Then the first-order expansion 
in position operator give; 



,15,25 



V = — F s • (rs z ) 



(7) 



with a spin force F s = — V (g* hbB(y)) applying on the 
carriers. Thus it becomes clear now that only when the 
spin current is defined as a time derivative of the spin 
displacement operator is z , can there naturally arises a 
conjugate driving force F s . As a consequence, the energy 
dissipation rate for the spin transport can be written as 
dQ/dt = J~ s -F s . It immediately suggests a thermody- 
namic way to determine the spin current by simultane- 
ously measuring the Zeeman field gradient (spin force) 
and the heat generation. 

The existence of a physical driving force for the con- 
served spin current makes it possible to construct On- 
sager relation between spin transport coefficients and 
other kinds of force-driving transport coefficients. For 
instance, since an electric force E may drive a spin-Hall 
current through the spin-orbit interaction, one naturally 
expects that a spin force F s may also induce a charge- 
Hall current. Naturally, an exact Onsager relation be- 
tween these intrinsic Hall effects may be established: In 
a general sense, when two kinds of different forces, say a 
spin force F s and an electric force E coexist as driving 
forces, then the linear charge-current and spin-current 
responses to these two forces can be expressed as 



Js 

J, 



a 
a 



a 



F s 
E 



(8) 



where J s is spin current and J c charge current, a ss 
(*ct* cc ) is the spin-spin (charge-charge) 3x3 conductivity 
tensor characterizing spin (charge) current response to a 
spin (charge) force F s (E). In the same manner, the off 
diagonal block < o >sc denotes spin current response to an 
electric field, and er cs denotes charge current response to 
a spin forced. A general relationship between a sc and 
a cs can be explicitly derived with a proper definition 
of the spin current, or imposed by the general Onsager 
relation 2 ^: 



T a0 



(9) 



where e a , £p are equal to +1 or —1 depending on whether 
the displacement (corresponding to current) operator is 
even or odd under time reversal operation T. In the 
present case, the displacement for a spin force is odd 
under T, while the displacement for an electric field is 
T-invariant, implying e a ep = — 1. 

Here the key point is that the Onsager relation is only 
attainable when the spin current is corrected by a torque 
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dipole term. To see this more clearly, we employ a stan- 
dard Kubo formula description as follows: Let us put a 
spin force along /i-direction and an electric field E along 
^-direction on equal footing by including both of them in 
the total Hamiltonian, H = Hq — Fid\ — F2CI2, where 
F\ = E and F2 = F s are the generalized forces ap- 
plied on the charge and spin degrees of freedom, whereas 
di = —ex^ and c?2 = s z x v are the corresponding displace- 
ment operators in which x^ denotes the /x-component of 
the position observable r. The response currents are ob- 
tained as the expectation values of the generalized veloc- 
ity operators in the perturbed states, 



di = -eiu, d 2 = 



SzXv T" SzXiy 



(10) 



where the standard symmetrization procedure in d 2 is 
implied. The presence of these two external forces 
will change an arbitrary stationary quantum state 
\a) of the original Hamiltonian into \a') — \a) + 
E/3#a,i=i,2 \P){0\Fidi\a)/(e a - ep + irj), where e Q is un- 
perturbed energy for \a) and 77 — ^ arises from the fact 
that when time t — ► —00, the perturbation is adiabat- 
ically switched off. Then the currents as a linear re- 
sponse to the perturbations are, after a straightforward 
derivation and by noting the substitution (f3\di\a) = 



ih 



given by (di) = J2 a ' /<*' 



y^,j cTjjFj with the conductivity matrix 

U V = ~ 1 „2 U« - IP) l 11 ) 



Q#/3 

E 



(e Q - e/3) 2 + V 2 



df a (e) 
de 



Re[(a\d 4 \P){l3\d j \a)}d(e a -ep) 



Here /„ is the equilibrium fermi function for band e a . 
Obviously, the first term is antisymmetric under the 
inter-exchange i <-> j, while the second term is the sym- 
metric part of Oij . While the symmetric part denotes dis- 
sipative contribution to the charge or spin transport, the 
first antisymmetric term is dissipationless and is relevant 
to our discussion in this paper. For an ideal system with- 
out scattering and many-body interaction, it is obvious 
that the dissipative (symmetric) part in Eq. (|lip vanishes. 
The dissipationless conductivity has manifested itself in 
a fundamental way in condensed matter physics involv- 
ing such issues as quantum and anomalous Hall effects. 
Our present discussions, as well as recent active studies 
of spin-Hall effect, are also focused on this dissipationless 
part of the transport conductivity. 

It becomes clear from Eas. pll)) - ([TT|) that the intrinsic 
spin-Hall conductivity is given by 



Hlm[(c 



aj M |/3></3|i„|a>] 



[e a - ep) 2 + rj 2 



where again symmetrization is implied in the product 
of two operators. On the other side, the dissipationless 



inverse spin-Hall conductivity driven by a spin force is 
given by 



E 



hIm[(a\x v \P)(p\s z x^ + s z x^\a)] 



(£ a - £p) 2 + V 2 



(fa ~ ff>) 
,( 13 ) 

The antisymmetric property shown above immediately 
gives the equality 



(14) 



which coincides with the general Onsager relation J9J. 
Note that Eq. (fT4]) denotes the Onsager relation only for 
the dissipationless (antisymmetric) part in conductivi- 
ties. A full Onsager relation by taking into account impu- 
rity scattering and an external uniform magnetic field Bo 
turns out to be given by a^ u (B ,T ) = -a^^-Bo, -to), 
where To characterizes scattering-induced relaxation 
time. Therefore, from Eqs. p^)l - (fT5|) one can see that 
to ensure the Onsager relation between the two kinds of 
Hall conductivities, the spin current should be defined as 



Js 



dt 



(15) 



instead of conventional definition. 

The above arguments clearly show two prominent fea- 
tures of the new spin current which is absent within the 
conventional definition of the spin current: (i) The spin 
current is now conserved with a physical conjugate driv- 
ing force; (ii) The Onsager relation is built up. Besides 
these two prominent featues, there is another physically 
valid property: For simple insulators whose single par- 
ticle eigenstates are localized (Anderson insulators), the 
spin transport coefficients vanish. Indeed, for spatially 
localized eigenstates, we can evaluate the (intrinsic) con- 
ductivity from Eq. (TT2]) as, 

= -e^/ Q , _ )2 (16) 

= -eh^2 f a (a\[rs z ,r]\a) = 

a 

where we have used (a\d(rs z )/dt\(3) — j L { e a — 
ep)(a\rs z \P) and </3|v|a) = £( C/ j - e a )(fi\i \a). The in- 
volved matrix elements are well defined between spatially 
localized eigenstates. 



IV. KUBO FORMULA FOR CONSERVED SPIN 
TRANSPORT COEFFICIENTS 

In this section we show how to evaluate in practice the 
conserved spin Hall conductivity based on the new defini- 
tion of spin current in a crystal. A formal description has 
been given in Eq. (10) for the conserved spin conductiv- 
ity, while the general states \a) should now be replaced 
by the electron Bloch wave function ip n k(r). Clearly, 
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the conserved spin Hall conductivity includes two com- 
ponents, 

a % = G fv + O^v ( 17 ) 

The first one is the usual conventional spin Hall coeffi- 
cient, which is ready to be rewritten in the Bloch-state 
space as 

C = -eft E t/( e « k ) ~ /( £ ™' k )] x ( 18 ) 
Im(u nk |^{'0 M ,s z }|u n / k )(u n / k |^|Mn k ) 

(e nk ~ £n'k) 2 + f/ 2 

Here v(k) dH(k)/hdk, H(k) 

exp(— ikr)H exp(ik-r), |u„ k ) is the periodic part 
of the electron carrier Bloch wave function, and (— e) 
is the electron charge. The Kubo formula (fl~8|) for the 
conventional spin Hall conductivity has been used by 
most of previous investigations^. 

The second component crf£ in the conserved spin-Hall 
conductivity, as shown in Eq. ^ and Eq. fT2")) , comes 
from the contribution of the spin torque density term. 
Due to the presence of the position operator in this term, 
and the fact that the position operator is not well-defined 
in the Bloch-state space, one needs to take some care in 
dealing with o s * v . As one choice for derivation, we pro- 
ceed by a regularization scheme, s^x^ — > s z sin^x^)/^, 
followed by the limit — ► 0. The symmetrized /i- 
component of the spin torque density operator occurred 
in Eq. (|12p can now be rewritten as 

Pr = ks z , S J^1} = lim ± [{s z ,e^} - {s z ,e-^}] 
2 q q^a 4iq 1 J 

(19) 

Substituting Eq. (fTO|) into Eq. (JT2J) and after a straight- 
forward manipulation, we obtain the Kubo formula for 

a% = -eh\im - £ [/(e nk ) - /(e„, k+q )] X (20) 
Re[(n nk |f(k,q)|u 7I > k+ q)(M n > k+ q|v(k,q)|u 7lk )] 

(fink - £n'k+q) + "T\ 2 

where f(k, q) = \ [f(k) + f(k + q)] with f (k) = 
(l/ih)[s z ,H(k)}, v(k,q) = i[v(k)+v(k + q)]. The limit 
of t] — > should be taken at the last step of calculation, 
and as a result, there is no intra-band (n — n') contribu- 
tion. Note that to properly calculate <j s Z, in practice, all 
the terms in Eq. (4) with the subscript k + q should be 
expaned at k to first order in q. 

Equation (j2"0|) can also be derived from an analysis of 
the spin torque dipole density T z (r) which can be de- 
termined unambiguously as a bulk property within the 
theoretical framework of linear response. Consider the 
torque response to an electric field at finite wave vector 
q, T z (q) — x(q)-E(q). Based on Eq. ^ which implies 



7^(q) = — iq-P T (q), we can uniquely determine the static 
response (i.e., q — > 0) of the spin torque dipole: 

P r (q) = Re{ i V q [ X (q)-E]} q=0 . (21) 

Here we have utilized the condition x(0) = 0, i.e., there 
is no bulk spin generation by the electric field. Thus 
the consequent spin-transport coefficient for the electric 
response of the spin torque density is given by 

<7£=Bfi[i0,„x*(q)]q=o. (22) 

Again within the standard Kubo-formula formalism, the 
spin torque response coefficient x(q) can be calculated to 
be 

X(q) = ieh t/( £ »k) - /(en'k+q)] x (23) 

Re[(u nk |f(k, q)|Un'k+q)(Mn'k+q|v(k, q)|Unk)] 
(frik - Cn'k+q) 2 + ?? 2 

which, combining Eq. (|22p , gives the same expression for 
a^l as that in Eq. (|20J) . 

Due to the inclusion of contribution from the spin 
torque dipole term, one may expect that in some spe- 
cial cases the transport properties of the conserved spin 
current are essentially different from that of conventional 
spin current. This will turn to be true by studying in 
the next section the intrinsic spin Hall conductivity in a 
2DHG system, which is modeled by a combined Luttingcr 
and spin-3/2 SIA Rashba spin-orbit interactions. For dis- 
ordered systems, the spin-Hall conductivity based on our 
new spin current has also been calculated in a recent 
work^l, and is found to depend explicitly on the scat- 
tering potentials for the two dimensional Rashba models 
with /c-linear or fc-cubic spin-orbit coupling. For one who 
instead prefers the conventional definition of the spin cur- 
rent, the above discussions are of course still helpful due 
to the fact that one cannot at last avoid tackling the 
calculation of the spin torque contribution to the total 
transport spin current and the spin accumulation, while 
our above derivation clearly indicates how to do that in 
practice. 

V. APPLICATION TO 2DHG SYSTEM 

As a practical application of the above general theory 
of the conserved spin current to the real physical sys- 
tems, in this section we focus our attention to the intrin- 
sic spin Hall effect in a 2DHG system, which has been 
experimentally^ , and theoretically^ investigated within 
the conventional spin-transport framework. Following 
Refi2£, the spin-orbit interaction in this system is mod- 
eled by combined Luttinger and spin-3/2 SIA Rashba 
terms. The resultant kp Hamiltonian reads (h is set to 
be unity): 

H = ( 7l + \l2)^ - ^(k-S) 2 + a(S x k)-2, (24) 
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where the confinement of the quantum well in the z direc- 
tion makes the momentum be quantized on this axis. The 
crucial difference between SIA Rashba term in 2DHG and 
S1A Rashba term in the 2DEG lies in the fact that S 
in 2DHG is spin-3/2 matrix, describing both the heavy 



(HH) and light (LH) holes. For the first heavy and 
light hole bands, the confinement in a quantum well of 
thickness a is approximated by the relation (k z ) = 0, 
(fc 2 ) w (ir/a) 2 . The energies for HH and LH are given by 



2L(k 2 + (k 2 z ))±\ak- 
2m 2 



a 2 k 2 ± Hlk{k 2 + (fc 2 )) + -^U 2 , 
m m 



(25) 



7i 
2m 



(fc 2 + (k 2 z )) ± ^ak + \/a 2 k 2 ± ^-Hk 2 + (fc 2 )) + ^d 2 , 




k (nm" 1 ) a (meV nm) 



FIG. 1: (a) Approximate band structure of the GaAs 2DHG; 
(b) heavy-hole fermi wave vectors as a function Rashba coef- 
ficient a for rih = 1.5 x 10 12 cm -2 . 



which is in agreement wit h 29 i 30 . Figure 1(a) gives a typi- 
cal band structure for GaAs (71 = 6.92, 72 = 2.1) with a 
r point gap of 40 meV and a Fermi momentum splitting 
of the hole band at Fermu momentum (0.2 nm -1 ) of 5 
meV, which requires a SIA splitting a ~ 50 meV-nm. In 
the recent experiment of spin-Hall effect 10 , this energy 
gap is of order AE = 40mcV, which corresponds to an 
a = 8.3 nm thick quantum well. During the following nu- 
merical calculation, these three material parameters (71, 
72, and a) will be fixed to be the values mentioned above, 
whereas the Rashba coefficient a and the hole density 
are treated as the tuning parameters. 



where and in the following k 2 = k 2 + fc 2 , d = 

y/k 4 + (fefp - k 2 (k 2 ). The HH and LH bands are split 
at the r point by A = 272 (fc 2 ) /m. Depending on the 
confinement scale a the Luttinger term is dominant for a 
not too small, while the SIA term becomes dominant for 
thin quantum wells. 

By expanding the above formulas for small k << (k z ) 
it is seen that the spin splitting of the HH bands is fc 3 
whereas the spin splitting of the LH bands is k^, 

e hh _ e hh = 3 a(a 2 - 4%r(ff)) fc3 6 

8 4(fc z 2 ) 2 

e lh _ e lh = 2ak + 0(k 3 ) 



First we consider the simple case of thin quantum 
wells. In this case the SIA Rashba term can be ne- 
glected (a — 0) and a full analytic calculation an be 
carried out. In the absence of SIA term, the (normalized) 
eigenstates for HH and LH bands are given by uf H = 
ai||,f> + a 2 fe^||,-|), u HH = ai|§,-§) + a 2 k 2 __ 
u\ H = ai|§,±) - a a fe|f.|§, -§), and U - H = a ilf>-5) - 
a2fc 2 ||,§), where ||,m z ) (m 2 = §,...—§) are the eigen- 
states of S z , k± — k x ± iky, a% = y~(d — \k 2 + (k 2 ))/2d, 

a,2 = yj3/8d(d — ifc 2 + (fc 2 )). Then using the standard 

Kubo formula (fT9| for the conventional spin Hall conduc- 
tivity, it is straightforward to obtain a s ® as follows 



47T 



HH 



JLH) 



3^ 
d 2 



(4 



72 d 2 



aia 2 fc 2 (a 2 + alk^^i + 72) 



kdk. 



(27) 



where Jhh and Jlh are fermi functions for HH and LH tribution from the torque dipole term turns out to be 
respectivel. The spin Hall conductivity due to the con- 
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FIG. 3: (Color online) (a) conventional spin Hall conductance 
component crJJ°, (b) spin torque dipole component a y T x , and 
(c) the total spin Hall conductance a yx (in unit of e/4n) as a 
function of Rashba coefficient a and hole density nu (contour 
plot). 



FIG. 2: (Color online) conserved spin Hall conductivity a yx 
(black line), and its two components a y % (blue line) and a yx 
(red line) as a function of the heavy hole fermi wave vector 
for a = and n h = 1.5 x 10 12 cm -2 . 



The first line in Eq. (|28[) i s obtained by expanding 
|«n'k+q), f(k,q), and e„'k+ q in Eq. ((201) at k to the 
first order in q while remaining other quantities to be 
their values at q = 0. Whereas the second line in Eq. 
([28]) with df /dk is obtained by a linear expansion of the 
fermi function with respect to q. The other terms occur- 
ring in Eq. (|20|l turn out to take no contribution to a yx 
for the present model. In the strong confinement limit 
(k 2 < (fc 2 )) and at zero temperature, Eqs. 1(27 ) - (|2g ]) are 
reduced to 



so 

'yx 



e 

47T 

e 

-in 



3(fc 



4 

HH 



^lh) 



(7i+72)(fc| 



HH 



^lh) 



3 7 i(fc 



6 

HH ' 



4 7 2 (k 2 z 



^lh) 



3(fc 4 



HH 



^lh) 



4 72 (fc 2 



(29) 



where kun and k^n are the fermi wavevectors for HH 
and LH respectively. 

Since for the experiment data, the LH bands are un- 
occupied (by the holes), so during numerical calculations 
throughout this paper we have set the fermi wavevector of 
the light holes to be zero, k^n — 0. Based on Eqs. (f27|) - 
([28| , figure 2 shows the conserved spin Hall conductivity 
a yx (black curve) as a function of the fermi wavevector 

knn (scaled by the confinement wave vector \J (A; 2 )) of 
the heavy holes. For clarification, we have also plotted 



in Fig. 2 the separate contributions from the conven- 
tional part a yx (blue curve) and spin torque dipole part 
a yx (red cuve); their sum gives a yx . With increasing the 
fermi wavevector knn, one can see from Fig. 2 that the 
two components a yx and a yx increase in amplitude but 
always differ in a sign. The resultant total spin Hall con- 
ductivity <j yx displays a non-monotonic behavior. The 
typical experiments 1 ^ are usually carried out in the re- 
gion of small value of fermi wavevector knn (compared 
to y/ (fc 2 )). In this region, it reveals in Fig. 2 that the 
tendency of the total SHC aligns with that of the conven- 
tional one with very little difference. Under experimental 
condition = 0.38 nm and knn — 0.31, the con- 

ventional term is obtained to a yx = —^f" (compared 
to the previous 2 ^ theoretical calculation result of — 
while the spin torque dipole term is given by a yx — ^j^- 
The total spin conductance is therefore a yx = - 
a tiny deviation from the value of conventional one a yx . 
We notice that the numerical simulation 10 related to the 
experimental setup gives a similar value of a yx , although 
the experiment was done in strong SIA Rashba spin-orbit 
coupling regime, while the present calculation with the 
result given in Fig. 2 is for a = 0. 

Now we turn on the SIA Rashba term which is essen- 
tial for understanding spin transport in strong spin-orbit 
coupling regime in experimental quantum well setup. In- 
clusion of cv-term in makes the analytic derivation of spin 
Hall conductance very tedious. Instead of doing that, we 
numerically obtain the eigenstates and eigen-energies of 
Hamiltonian (j2"4")) , and substitute them in Eq. (|18p and 



with 

4-7T ' 

SO 
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(|20| . After a k- integration with fermi wavevector pinned 
by the given hole density at zero temperature [see Fig. 
1(b) for the relationship between HH fermi wavevectors 
and a], the spin Hall conductance is therefore steadily 
calculated in a wide range of system parameters. Figures 
3(a)-(c) summarize the behavior of the conventional, the 
spin-torque-dipole contributed, and the total spin Hall 
conductance, respectively, as a function of Rashba coef- 
ficient a and hole density nu (contour plot). Compared 
to the results without consideration of Rashba spin split- 
ting (Fig. 2), the prominent new features occurred in 
Figs. 3 are: (i) The conventional term a yx and the spin 
torque dipole term a yx jump to much larger values upon 
switch on of Rashba spin splitting. This jump is caused 
by the presence of coupling between the two HH bands, 
which is absent without considering the SIA term; (ii) 
The amplitude of <r yx becomes larger than that of con- 
ventional one a yx with a sign difference. In fact, the 
typical value of a yx in a wide range in (a, rih) parameter 
space is about ^^f-, while a yx is typically of re- 
sulting in dy-r^ff . Thus overall the conserved spin Hall 
conductance changes a sign compared to the previous 
calculations based on conventional definition; (iii) The 
stripe in Figs. 3 reveals that the amplitude of intrinsic 
spin Hall conductance for 2DHG is insensitive to the am- 
plitude of Rashba coefficient a, which implies a universal 
character. The previous studies^ on a fc-linear Rashba 
2DEG system within the conventional definition of spin 
current have revealed a universal spin Hall conductance 
(a yx ) of e/87r. We have also calculated the conserved spin 
Hall conductance <jf„ with the same model and found it 

y 

to be — e/87r. Together with the above results for 2DHG, 
one can see the key role played by the spin torque dipole 
term, which tends to overwhelm the conventional spin 
conductance by the opposite contribution. 

Since the LH bands are usually unoccupied by the holes 
in the experiments, and one can expect that in the weak 
Rshba spin splitting, the major contribution to the spin 
Hall conductance comes from the coupling between the 
two HH bands. In this case we can approximate the four- 
band Hamiltonian by an effective two-band one. Af- 
ter a standard adiabatic elimination procedure, we obtain 
an effective two-band Hamiltonian H for heavy holes (in 
the strong confinement limit k 2 < (k 2 )), 



H 



k 2 



(k 2 z ) 



ZiriHH 2m* 



i(3(k 3 _a^ 



kla-) 



(30) 



where /3 = ^ a({k *\ k t^ LHa \ m HHiLH = m/(7i T 272), 

to* = 2mHH'm<LH / (mim + ttilh) is the reduced mass 
which is due to the proper account for the coupling be- 
tween HH and LH bands, and a± — a x ±ia y are Pauli ma- 
trices. To reflect the angular momentum quantum num- 
bers of the heavy holes, the conserved spin Hall current 
operator is defined by J s = | d . Thus after an adia- 
batic projection onto the two HH bands, the Hamiltonian 
is reduced to an effectice fc-cubic Rashba model. In 
this case, we can derive an analytic expression for a yx 



and its two components a yx and <r yx . The results at zero 
temperature are given by 



s0 
'yx 



' 32tt to/3 [ k+ kj 



yx 32ttto/3V >-> °- 



27e 7x , 1 1 

k j 

9e 71 , 1 1 

167r to/3 fcy" k~f 



9e 
87 
9e 
8^ 



(31) 



where k~J and kj are fermi wavevectors of the two HH 
bands. According to the experimentally accessible hole 
density rih, these two fermi wave vectors can be expressed 



k + '~ - 



J_ 

,2/3 



1-1 1- 



2/3 s 



47m/! + 3irnh 



1/2 



1 2. 

2 2/3 



2/3 
7 



where 7 = h 2 /2m*. In the limit n h « (l/47r)(7/2/3) 2 , 
these fermi wavevectors have more simple forms: fct' ~ 
\/2imh T (2/3/7)7^. In this case one gets a yx = — |f , 

°yl = If' thus a U = If- In § eneral Vyx depends on 
the spin-orbit coupling and the fermi energy which is 
pinned by the hole density. For further illustration, we 
show in Fig. 4 the numerical results of conseved spin 
Hall conductance crfL as a function of Rashba coefficient 
a for the four-band (black curve) and the approximate 
two-band (red curve) Hamiltonians. Note that the coef- 
ficient (3 in Eq. (|30p is obtained from a by the equation 
below Eq. (f3"0"|) . It reveals in Fig. 4 that the results 
from the two models agree reasonably in the experimen- 
tally relevant variation of spin-orbit coupling coefficient. 
When the value of spin splitting akt (kf characterizes 
fermi wavevector) is large enough to be comparable with 
the confinement characteristic energy (k 2 )/2m, then the 
deviation of the two-band approximation from the four- 
band calculation becomes obvious, which is featured in 
Fig. 4 by the fact that while the two-band a yx increases 
monotonically as a function of a 2 (which can also be seen 
from Eq. (pTT]) , the four-band spin Hall conductance dis- 
plays a much weaker dependence of Rashba coefficient. 
Therefore we arrive at a conclusion that the neglect of 
HH-LH coupling is no longer valid for strong spin split- 
ting and a more exact multi-band simulation is necessary 
for the general 2DHG systems. 



VI. SUMMARY 

In this article we have studied spin transport features 
in general spin-orbit coupled systems. Due to the spin- 
orbit coupling, the spin is not conserved by the occur- 
rence of spin torque during its flow through the sample. 
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FIG. 4: (Color online) conserved spin Hall conductivity at x 
as a function of Rashba coefficients a for nj, = 1.5 x 10 12 
cm -2 . The black line is the exact four-band result while the 
red line gives the approximate two-band result for the heavy 
holes. 

Therefore, how to describe the spin transport within an 
intuitive drivc-thcn-flow picture becomes an important 
issue. The key point clarified in this paper is as follows: 
(i) We have shown within the linear response formalism 
the necessity to include the spin torque dipole term in 
the expression of the spin current. The real advantage 
of our new definition of spin current lies in the fact that 
it provides a satisfactory description of spin transport 
in the bulk. With our new spin current, one can now 
use the spin continuity equation ^ to discuss spin accu- 
mulation in the bulk, e.g., by generating a non- uniform 
electric field or spatially modulating the spin Hall con- 
ductivity. Our new spin current vanishes in Anderson in- 
sulators either in equilibrium or in a weak electric field, 
which enables us to predict zero spin accumulation in 
such systems. More importantly, it possesses a conju- 
gate force (spin force), so that spin transport can be 
fitted into the standard formalism of near equilibrium 
transport. The conventional spin current does not have 
a conjugate force, so it makes no sense even to talk about 
energy dissipation from that current. The existence of a 
conjugate force is crucial for the establishment of On- 
sager relations between spin transport and other trans- 
port phenomena, and its measurement will be important 
to thermodynamic and electric determination of the spin 
current, (ii) A general Kubo formula for the conserved 
spin transport coefficients, consisting of the conventional 



and spin-torque-dipole contrbutions, has been derived, 
which makes the practical bulk calculation to be feasi- 
ble. 

Based on the Kubo formula for the conserved spin cur- 
rent, we have analyzed in detail the spin Hall effect in 
2DHG system with the parameters chosen to be rele- 
vant to recent experimental measurement on GaAs quan- 
tum well and modeled by the Hamiltonian consisting of 
a Luttinger spin-orbit coupling term and a S1A Rashba 
term. In the absence of Rashba spin splitting, the two 
HH (and LH) bands are degenerate and the non-zero spin 
Hall conductance only comes from HH-LH transition. In 
this case, we have derived an analytic expression for the 
conserved spin Hall conductance a yx and its two com- 
ponents, the conventional one a yx and the spin torque 
dipole correction a yx . These two components have been 
shown to compete each other, which is verified by a dif- 
ference of sign between them. In the case that only the 
Luttinger term is taken into account, it has been found 
that the amplitude of a yx is much small than that of a yx . 
In fact, the value of a yx is typically of — while a yx is 
about Sp 2 . Thus in this case, the spin torque dipole cor- 
rection to the spin Hall conductivity is relatively small. 
When the SIA Rashba spin splitting is turned on, we have 
found that there occurs a jump in amplitude for both a s yx 
and a s yx . a yx changes now to take a typical value of — 
while a yx jumps to a characteristic value of -jjp There- 
fore, the presence of Rashba term not only stirs up a large 
enhancement of the total spin Hall conductivity and its 
two components, but also changes a fundamental sign for 
a yx because the spin torque dipole correction now over- 
whelms the conventional contribution. This jump and a 
change of sign for a s yx comes from the HH-HH coupling, 
and therefore can be modeled by an effective two-band 
Hamiltonian which has been done in this paper by adia- 
batically eliminating the LH states and reducing the sys- 
tem to a fc-cubic Rashba model with a dressed spin-orbit 
coupling coefficient {(3 in Eq. ([30]) ]. By a detailed com- 
parison, furthermore, we have shown that this two-band 
a yx displays a quadric dependence of original Rashba co- 
efficient a, while the exact four-band treatment shows a 
somewhat universal character for a yx . Thus the neglect 
of HH-LH coupling is no longer valid for strong spin split- 
ting and a more exact multi-band simulation is necessary 
for the general 2DHG systems. 
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